function F = LorentzForceCase1_Func_Z(theta,r,l)
X = r.*cos(theta);
Y = r.*sin(theta);
Z = l;
I = 1;
dw = 0.35e-3; 
Area = pi*(dw^2)/4;
J = I/Area; 
ao_em = 30e-3; 
ao_pm = 30e-3;
ai   = 13e-3;
l_em    = 10e-3;
delta_h = -10e-3;
l_pm = 6e-3;
M = 1.465;
B3d_EM_Cir = arrayfun(@(X,Y,Z)J*MagneticField3D_Cir(X,Y,Z,l_em,ai,ao_em));
B3d_PM_Cir = arrayfun(@(X,Y,Z)CylPM_flux_Cir(M,ao_pm,l_pm,X,Y,delta_h+Z));
B3d = B3d_EM_Cir+B3d_PM_Cir;
B3d_x = B3d(1);
B3d_y = B3d(2);
 F = (r.*(-sin(theta)).*B3d_y - r.*cos(theta).*B3d_x;
for 

end

